Quasilinear theory of the 2D Euler equation 



Pierre-Henri Chavanis 

Lahoratoire de Physique Quantique, Universite Paul Sabatier, 118 route de Narhonne 31062 Toulouse, France 

(February 2, 2008) 



o ■ 
o ; 
o 

<: 



I , 
^ . 

q: 

O ■ 

^ ■ 



> 
O 



a^ : 

OS . 
o ■ 

• I— I 
>^ 

^ : 

Oh. 



X 



We develop a quasilinear theory of the 2D Euler equation 
and derive an integro-differential equation for the evolution of 
the coarse-grained vorticity Zj(r, t). This equation respects all 
the invariance properties of the Euler equation and conserves 
angular momentum in a circular domain and linear impulse 
in a channel. We show under which hypothesis we can de- 
rive a H-theorem for the Fermi-Dirac entropy and make the 
connection with statistical theories of 2D turbulence. 
PACS numbers: 47, 47.10.-)-g, 47.27. Jv, 47.32.Cc 
Two-dimensional flows with high Reynolds numbers 
have the striking property of organizing spontaneously 
into large-scale coherent vortices Q. The robustness of 
Jupiter's Great Red Spot, a huge vortex persisting for 
more than three centuries in a turbulent shear between 
two zonal jets, is probably related to this general phe- 
nomenon. Many other vortices are observed in geophysi- 
cal and astrophysical flows. Understanding the structure 
and formation of these organized states is still a challeng- 
ing problem. 

To be explicit, we consider as an initial condition a 
stripe of uniform vorticity to — ao surrounded by irrota- 
tional flow oj — 0. This stripe is unstable and generates a 
complicated mixing process leading to the formation of a 
quasi-stationary vortex slightly diffusing with viscosity. 
This is the classical shear layer (or Kelvin-Helmholtz) 
instability investigated numerically in, e.g., These 
authors propose to interprete the quasi-equilibrium state 
as a state of maximum entropy under the constraint of 
a fixed energy and circulation. This is motivated by the 
statistical theory of the 2D Euler equation developed by 
g Miller (1990) and [§ Robert & Sommeria (1991). A 
coarse-graining procedure is introduced and a mixing en- 
tropy is constructed to describe the chaotic interchange 
of vorticity levels along the evolution. Since the vorticity 
levels cannot overlap, they follow an exclusion principle 
and this leads to a statistics of a Fermi-Dirac type. Com- 
parision with numerical simulations shows a very good 
agreement with the theoretical prediction in the core of 
the vortex where the fluctuations are sufficient to validate 
the ergodicity hypothesis. This statistical mechanics of 
phase mixing is closely related to the theory of "violent 
relaxation" developed by Lynden-Bell (1967) for colli- 
sionless stellar systems (e.g., elliptical galaxies) described 
by the Vlasov equation 

Less is known concerning the relaxation towards equi- 
librium. This is clearly a complicated task and analytical 
results will be obtained only by introducing approxima- 
tions. Our objective is to derive a kinetic equation re- 



specting all the conservation laws and invariance proper- 
ties of the Euler equation and driving the system towards 
the Fermi-Dirac state by increasing the mixing entropy. 
If such a program can be realized this will provide a use- 
ful subgrid scale model allowing Large Eddy Simulations 
(L.E.S.) of 2D turbulence with potential applications in 
geophysics |^ . A first step in this direction was made by 
Q Robert & Sommeria (1992) using a variational proce- 
dure. They assumed that "out of equilibrium, the system 
evolves so as to maximize the rate of entropy produc- 
tion S while respecting all the constraints of the Euler 
equation" . This Maximum Entropy Production Princi- 
ple (MEPP) leads to an equation for the coarse-grained 
vorticity of a generalized Fokker-Planck type which can 
be compared succesfuUy with Direct Navier Stokes sim- 
ulations §,|l^. Their method was extended by jl^ Cha- 
vanis & Sommeria (1997) who derived a set of equations 
respecting, in addition, the invariance properties of the 
Euler equation. However, the MEPP is relatively ad hoc 
and assumes that the system evolves towards a maximum 
entropy state. In this letter, we consider a completely dif- 
ferent approach based on a perturbative expansion of the 
Euler equation. This is the counterpart of the quasilin- 
ear theory introduced in plasma physics and in stellar 
dynamics for the Vlasov equation [p"2|-p^. Within some 
approximations, we derive a new kinetic equation of a 
generalized Landau type for the coarse-grained vorticity 
and prove a H-theorem for the Fermi-Dirac entropy (in- 
stead of postulating it). The results of the MEPP are 
recovered as an approximation of our model. 

For a two-dimensional incompressible and inviscid 
flow, the Euler equation can be written: 



— + uVw = 
at 



u = -z A VV' 



u) = -Alp 



(1) 



(2) 



where luz = V A u is the vorticity and ■0 the streamfunc- 
tion (z is a unit vector normal to the flow) . The velocity 
can be expressed as an integral over the vorticity field as 



u(r,i) = J (fr'V{r' r)w(r',t) 



(3) 



where 



V(r'-.r).-i-(^+V.(r'^r) (4) 

represents the velocity created in r by a vortex of unit 
circulation located in r' (r^ is the vector r rotated by 



1 



+ ^). The term Vh(r' r) accounts for boundary effects 
(Vfo = in an infinite domain) and can be calculated with 
the method of "images" . 

In the situation described previously, the Euler equa- 
tion builds up an intricate filamentation at smaller and 
smaller scales. If we subdivise our domain into a lattice 
of macrocells of size e, only the "coarse-grained" vortic- 
ity uJ(r, t) can reach a stationary state. This "coarse- 
grained" vorticity is defined by a double averaging pro- 
cess [0| : a space average over the cell of size centered 
on r^'^and a statistical average to express our ignorance 
of the precise manner in which the phase filaments of 
vorticity ctq are distributed in the macrocell: 



(5) 



The fluctuating vorticity a) = cj — ZU, satisfying oj — 0, 
is simply the difference between the exact vorticity and 
the smoothed-out vorticity. The passage from discrete to 
continuous variables for to requires a hypothesis of scale 
separation. We shall assume that the velocity field u(r, t) 
consists of a strong large-scale component and a weak 
small-scale component such that the characteristic scale 
of u is much greater than that of u (we also assume that 
Lj ^ U). This hypothesis of scale separation was previ- 
ously made by ||l5| DubruUe & Nazarenko (1997). This is 
of course an idealization since the energy spectrum never 
presents a clear-cut gap in practice. However, this ap- 
proximation should reasonably well account for the non- 
local interactions between large eddies and small scales 
fluctuations in 2D turbulence. A more general study re- 
laxing this hypothesis can be found in fl^ ]. 

If we take the local average of the Euler equation (|^), 
we obtain a convection-diffusion equation: 



'dt 



uVZU = -VJ 



(6) 



for the coarse-grained field with a current J ~ uju related 
to the correlations of the fine-grained fluctuations. In 
turn, the fluctuations depend on the smoothed-out field 
according to the equation: 



— — h uVuj = — uVcj — UVUJ - 
dt 



uSJCj 



(7) 



obtained by substracting (|^) and (^. Within the scale 
separation hypothesis, we can neglect the non linear 
terms uVo) and uVo; which represent the interactions of 
the small turbulent scales among themselves How- 
ever, unlike [T^ , we keep the linear term uVtU which 
takes into account the interactions between small and 
large scales. Its order of magnitude uPe/L (where L ^ e 
is the domain size) is relatively small but this term has a 
cumulative effect (see equation |l^) giving rise to a diffu- 
sion process. We consider therefore the coupled system 



-;- h LUO = — VCJU 

at 



duj 
'dt 



Llu 



-uVw 



(8) 



(9) 



where L = uV is an advection operator. This "quasilin- 
ear approximation" is standard in plasma physics and in 
stellar dynamics for the Vlasov-Poisson system p^-[T^] 
but, to our knowledge, it has never been applied to the 
2D Euler system. Owing to the various approximations 
introduced, this theory can describe only the late quies- 
cent stages of the relaxation when the fluctuations have 
weaken. 

Introducing the Greenian: 



G{t2,ti 



exp< 



dtL{t) 



(10) 



we can immediately write down a formal solution of (|y), 
namely: 

Co(r,t) = G(t,0)(I;(r,0) 
- / dsG{t,t- s)u{r,t- s)\'Lj{r,t- s) (11) 

Although very compact, this formal expression is in fact 
extremely complicated. Indeed, all the difficulty is en- 
capsulated in the Greenian G{t, t — s) which supposes 
that we can solve the smoothed out Lagrangien flow: 



dr 



(12) 



between t and i — s. 

The objective now is to substitute the formal result 
( |TT| ) back into (js]) and make some closure approximation 
in order to obtain a self-consistant equation for ZZ)(r,t). 
If the vorticity were purely adverted by the stochastic 
velocity field u (like a passive scalar), the interaction (|^) 
would be switched off and we would end up with a diffu- 
sion equation for ZJ with a diffusion coefficient D ~ jTU^ 
where r is the decorrelation time |H,|l^. However, in the 
case of the Euler equation, the velocity fluctuations are 
induced by the fluctuations of the vorticity itself accord- 
ing to: 



u{r,t) = A / d^r'Y{r' r)w(r',t) 



(13) 



Therefore, considering (gT|) and (|13D, we see that the vor- 
ticity fluctuations a)(r,t) are given by an iterative pro- 
cess: uj{t) depends on u{t — s) which itself depends on 
uj{t — s) etc... Since |u|, of order ZUe, is much smaller 
than |u|, of order Luj, we can solve this problem pertur- 
batively. This is the equivalent of the "weak coupling 
approximation" in plasma physics [p^-p^. For conve- 
nience, we have introduced a counting parameter A in 
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( p^ which will be set equal to one ultimately. To order 
A^, we obtain after some rearrangements: 



xG'{t,t- s)G{t,t~ s) 



x{ Vir" r)w(r',i- s)w(r",i- s) — (r,i- s) 

duJ I 



+y'^(r" ^ r')c^(r,t - s)(i(r",t - s) — (r',t- s)| (14) 

In this expression, the Grecnian G refers to the fluid par- 
ticle r(t) and the Greenian G" to the fluid particle r'(i). 
The contribution proportional to A (not written) can be 
calculated with the assuption that uj is purely advected 
by the large scale velocity, i.e uVZU is neglected in (||). 
This is the case considered by However, in this ap- 
proximation the coarse-grained enstrophy J oJ^cPr is con- 
served [|l^ and no trend towards a self-organized state 
(e.g. maximum entropy or minimum enstrophy state) is 
apparent. The exchange of enstrophy between small and 
large scales (and also the source of entropy) corresponds 
to higher order corrections in the equation for uj. In this 
letter, we shall consider exclusively the term of order A^ 
which accounts for a diffusion process but we do not claim 
that the term of order A must be necessarily discarded. 

To close the system, it remains to evaluate the corre- 
lation function uj(r,t)oj{r' ,t). We shall assume that the 
scale of the kinematic correlations is small with respect 
to the coarse-graining mesh size and take: 



c;)(r, Oc;)(r', t) = e^6{r - r')i:}^{r, t) 



(15) 



This assumption is consistent with our scale separation 
hypothesis and was made previously by for the Eu- 

ler equation and by p2-n4| in plasma physics. Now: 



a)2 = (w — uj)'^ = Lip- — up' 



(16) 



For the case that we consider, the exact vorticity field uj 
can take only two values w = and w = uo • This implies 
that up — 10 — ctquj and therefore: 



(D(r, t)uj{r', t) — e^5{r — r')uj{aQ — uj) 



(17) 



Substituting this expression in equation (|IJ) and carry- 
ing out the integration on r", we obtain: 

duj 



xiVir'-^rP'iao-uj') — 



+V{r ^r')uJ{ao-Ij) 



dij' 
dr'" 



(18) 



We have written = uj{r' (t — s),t — s), Ut-s =uj{r{t- 
s),t- s), F'^(r' ^ r)t = Vir'it) r{t)) and V{r' - 



r)t-s = Vir'it - s) ^ r{t - s)) where r{t - s) is the 
position at time t — s of the fluid particle located in r = 
r(t) at time t. It is determined by the characteristics ( p^ ) 
of the smoothed-out Lagrangian flow. 

Equation (^8|) is a non Markovian integro-differential 
equation: the value of tU in r at time t depends on the 
value of the whole vorticity field at earlier times. If the 
decorrelation time r is short, we can make a Markov 
approximation and simplify the foregoing expression in 



duJ e^T d 

— h UVUJ — 

dt 2 drf' 



J d2r'y'^(r' ->r) 
x|y''(r' ->r)tZJ'(ao -U') 
+V''{r^r')uJ{(Jo-Zj) 



duJ' 



dr 



In the case of an infinite domain, V(r r') 
r) and we have the further simplification 



(19) 



-V(r' 



duj 
'dt 



87r2 drf^ 

) — 
dr 



-I I -/N duj _ duj' 
^('^o-w)7^-^(^o-^)^ 



where 



K^"'{v' - r) = 



(20) 



(21) 



and ^ = r' — r. The symmetrical form of this equation 
is of course reminiscent of the generalized Landau equa- 
tion in plasma physics obtained with a quasilinear theory 

Introducing a tensor 



and a vector 

^2 



d^'V^ir' ^ r)V"{r' ^ r)w'(ao - ZJ') 

(22) 



r,'' J d'^r'V'ir' r)V''{r r') 



duj' 



dr'' 



(23) 



equation ( [l9| ) can be rewritten in the more illuminating 
form: 



duj d 



duj 



D^^'^—+u;{ao-ZJ)v^ 



(24) 



This equation has the structure of a generalized Fokker- 
Planck equation with a diffusion term and a drift term. 
The diffusion term corresponds to the turbulent viscos- 
ity introduced ad hoc in most parametrizations of tur- 
bulence. However, this term alone breaks the conser- 
vation laws of the Euler equation. The present theory 
shows that an additional drift term must exist in order 
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to recover these properties. The drift is non hnear in 
uj so that ( p^ ) is not, strictly speaking, a Fokker-Planck 
equation. This nonhnearity accounts for the constraint 
(^(r, t) < ctq imposed at any time by the conservation of 
the fine-grained vorticity (see equation (|l|)). 

Equation ([l^) respects the invariance properties of the 
2D Euler equation and has the same structure as equa- 
tion (23) of [0 Chavanis & Sommeria (1997) derived 
on the basis of thermodynamical arguments. In their 
work, the constraints of the Euler equation were satisfied 
with the aid of Lagrange multipliers. In this new ap- 
proach, the conservation laws follow naturally from the 
symmetrical structure of equation (|l9|), as for the usual 
Landau equation (the linear impulse P = J ujydP'v and 
the angular mometum L = J ujr'^cPr play the role of the 

impulse P = / /vd^v and kinetic energy K = J f^d^-v 
in plasma physics). This is more satisfying from a phys- 
ical point of view. Moreover, in the thermodynamical 
approach, the increase of entropy is postulated whereas 
in the present situation a H-theorem for the Fermi-Dirac 
entropy 



S=- 



l^ln^ 



1 llnf 1 



d^r (25) 



^0, 

results immediately from equation (|19| ) . This is proved 
by taking the time derivative of (|25|), substituting for 
(p^), interchanging the dummy variables r and r' and 
summing the two resulting expressions. Of course, the 
increase of entropy is due to the coarse-graining proce- 
dure which creates some irreversibility (the indetermina- 
tion on the position of the vorticity levels in a cell). The 
entropy S'[w] for the exact vorticity uj is conserved by the 
Euler equation as the integral of any function of uj. 

It is remarkable that a quasilinear theory is sufficient 
to generate a turbulent viscosity (but also a drift) and a 
source of entropy. We do not necessarily have to advocate 
the non linear terms in (|^) to get these properties. Note 
also that the entropy associated with the (coarse-grained) 
Euler equation is the Fermi-Dirac entropy (|2^) in agree- 
ment with the works of at equilibrium. Unfortu- 
nately, equation (^|) does not conserve energy exactly. 
Therefore, the system will ultimately relax towards the 
solution ZZJ = ao / {1+ Xe'xp{aaQr'^)) which is the maximum 
entropy state at fixed circulation and angular momen- 
tum. This means that our approximations break down 
at very late times. 

A further connection with the statistical theory of 2D 
turbulence can be found. Equation (|l^) is an integro- 
differential equation whereas the equations derived from 
the MEPP ||,|,|ll| are diflFerential equations. The usual 
way to transform an integro-differential equation into a 
differential equation is to make a guess for the function U' 
appearing in the integral. It makes sense to replace ZU' by 
its optimal value uj' — 170/(1 + '^exp(/?croV'')) maximizing 
entropy at fixed energy and circulation. Substituting in 
(E2) (E3) and making a "local approximation" , we obtain 



ri = D(3ViIj 



(26) 



(27) 



In equation ( |2^ ) , we recover the form of the drift derived 
by |l| Chavanis (1998b) in a point vortex model. The 
drift coefficient can be interpreted as an Einstein formula. 
Substituting for the drift in (M) we recover the equation 



duj 

— -f uVuJ 

ot 



V{D{Vuj + l3ujiao - w) V^)) (28) 



obtained by Robert & Sommeria (1992) using a Max- 
imum Entropy Production Principle. Equation ( p8|) can 
be interpreted as a generalized Fokker-Planck equation 
|l7| . Note that the present approach provides the value 
( P7| ) of the diffusion coefficient which was left unknown 
by the variational principle |^ . This value coincides with 
the estimate of [^,|l^ based on a passive scalar model. 

The results of this letter can be extended to an arbi- 
trary spectrum of vorticity levels for a wider class of 
initial conditions. These results also complete the anal- 
ogy between 2D turbulence and stellar systems investi- 
gated by the author 
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